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Abstract 



Exploding wires are widely used in many experimental setups and pulsed power 
systems. However, many aspects of the process of wire fragmentation still remain 
' r~| '' unclear. If the current density is not too high, the wire may break up in the solid 

state. The experiments have shown that the wires break in tension due to longitudinal 
c/2 | forces of unknown nature. 

In a series of papers, see P, ^, ||, || Graneau argued that neither mechanical vi- 
brations induced by the electromagnetic pinch force nor thermal expansion could have 
^ | been responsible for the wire disintegration because they were too weak. To explain 

the phenomenon he appealed to the obsolete Ampere force law as opposed to the con- 
ventional Biot-Savart force law. Graneau argued that the Ampere force law would 
lead to a longitudinal tension in the wire, although his calculations may have been 
Q_i' in error on this point. Therefore, Graneau's explanations induced a controversy in 

electrodynamics with a number of authors arguing pro and con. 

Previous theoretical and numerical investigations served to provide a search for 
these forces without recourse to unconventional electrodynamics have identified the 
^vq ■ pinch effect and thermal expansion as a source of strong longitudinal vibrations. But 

<^> . the tensile stress component has been proved to be present only in the wires having free 

ends. Thus the mechanism doesn't give a satisfactory explanation of the phenomenon 
■ in the wires with clumped ends. In this investigation, we use a simplified magneto- 

thermo-elastic model to study flexural vibrations induced by high pulsed currents in 
c/2 \ wires with clumped ends on account of their role in the disintegration process. Several 

• ' aspects are studied, namely (i) the buckling instability due to simultaneous action of 

the thermal expansion and the magnetic force, (ii) the flexural vibrations induced in 
', initially bent wires. It is shown that the induced flexural vibrations are strong enough 

to lead to the breaking of the wire in a wide range of parameters. 



1 Introduction 

The phenomenon of wire fragmentation in the solid state by high pulsed currents was studied 
experimentally by Nasilowski [|IJ, and Graneau 0, || |J. They observed that a sufficiently 
strong electric current would shatter a thin metal wire under a broad range of conditions 
(various wire material and geometry, different current types, etc.). As a result of the explo- 
sion the wires fragmented into 2-100 pieces with apparent signs of longitudinal tensile stress. 
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The experimental results, as well as (sometimes controversial) early attempts to explain the 
phenomenon have been reviewed by Graneau || , Hong || and Molokov&Allen 0. 

It is obvious that an electric current passing through a metal wire induces stress waves. 
Their origins are in (i) the thermal expansion owing to the volumetric Joule heating, and (ii) 
the Lorentz-, or pinch-, force. These forces depend on the magnitude and the rate of change 
of the current passing through the wire. They grow as either of these parameters increases. 

Ternan || suggested that in Graneau's experiments with wires with free ends standing 
stress waves may be induced as a result of thermal expansion. Using a simple one- dimensional 
model he showed that the resulting stresses were sufficient to lead to a fracture. 

This mechanism of fragmentation of wires with free ends has been explored in more detail 
by Molokov&Allen in . They employed a magneto-thermo-elastic model of the stress- wave 
propagation, and solved the problem numerically on assumption of axisymmetric nature 
of vibrations . Their conclusion was that the characteristic value of both compressive and 
tensile longitudinal stresses obtained for the aluminium wire with radius a = 0.6 mm carrying 
current / = 5kA and having free ends (the parameters relevant to the experiment in ||) 
would be about 33 MPa per unit length of the wire. The stress has been shown to grow 
linearly with the wire length and thus, for sufficiently long wires, the resultant stress exceeds 
the ultimate strength of aluminium. This is 75 MPa at 100° C going down to 17 MPa at 
320° C §. 

Although the mechanism does give large value of tensile stresses, it cannot explain the 
fragmentation of wires with clamped ends (as in the experiments |l|, |]]) because within the 
axisymmetric model in this case the stresses can be only compressive J7|. It is clear, however, 
that in the wires with clamped ends flexural vibrations may be induced. The apparent signs 
of these modes have been observed in the Graneau's experiments on the wires with firmly 
clamped ends. The resultant wire shape, after breaking, clearly indicated on that, see 
photographs in |4[]. 

In this paper we carry out a special study of flexural vibrations to account for their role 
in fragmentation of wires with clamped ends. In section 2, we present the formulation of 
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the problem and discuss simplifications of the model. In section 3, the linear stability of the 
system is considered. In section 4, we carry out a qualitative analysis of expected stresses. 
In section 5, numerical algorithm and results are presented. 



2 Formulation 



Consider a wire which is firmly clamped between two points located along the z-axis. Let the 
distance between these points be L. We restrict our attention to the case of plane motion, 
and denote the deflection of the wire from the axis by X which is a function of z and time 
t. The deflection of the wire in the initial undeformed state is denoted by X . Possible 
three-dimensional modes, not considered here, may only increase stresses, especially at the 
clamped ends. 

Our theoretical consideration is based on a simplified model valid for sufficiently long or 
sufficiently thin wires. In this model the deflection X of the wire is assumed to be much 
lower than the wire length, i.e. \X\/L <C 1. Since large deflections of the order of the wire 
length are not expected in the wire explosion experiments, this restriction is not expected to 
be very stringent. Another small parameter used in the model is the ratio a/L « 1, where 
a is the wire radius. This condition is usually well satisfied. Plastic deformations are not 
considered here although they may be important in the vicinity of the melting point. Thus 
all deformations are assumed to be elastic. 

Then the deflection X obeys the following equation ( [|1C[| ) : 



n^ x + f q2 f x f 92 T d2x ° r d2x f dG * - n m 

where p is the material density, S is the wire cross-section which may be a function of z, 
J is the effective moment of inertia of the wire cross-section (for a circular cross-section 
J = 7ra 4 /4, where a is the wire radius), E is the Young modulus, G is the force along the z 
axis applied at the wire ends, F x is the distributed external force per unit length along the 
x direction and C y is the distributed moment of force along the y direction per unit length. 
For the physical situation considered here, the force G can be represented as 



G = ES{{1- l )/L - a[T(t) - T }}, 

The second term is a compressive force due to the thermal expansion, a being the linear ex- 
pansion coefficient. The first, nonlinear term is a retarding force due to the wire deformation 
from the initial state, where 

1 f L (dX\ 2 (dX^ 2 



is the increment of the wire length due to the deflection. When X = Xq, I = Iq and the 
retarding force equals to zero, as required. 

The boundary conditions are those for the clamped ends, namely 

dX 

X = — — = at z = and z = L. (3) 

oz 

The initial conditions are 

X = X t = 0. (4) 

2.1 The magnetic force 

In long, thin wires under consideration a/L << 1, uniform current is directed perpendicular 
to the cross-section. Thus any possible forces of magnetic nature act only in the cross-section 
plane. These forces create no moment and thus C y = in our case. 
The net distributed force F can be calculated directly from 



JjxBdS. (5) 



The integral in (|]) is taken over the wire cross-section area. The magnetic field is given by 
the Biot-Savart law 



* = T L !'^dV\ (6) 
Atc J R 3 



where dV' is the volume element on the wire, j is the current density, and the vector R = r— r 
is directed to the field point as usual. 

To calculate the x component of the force, F x , one can split the integral (^) over the wire 
length z into two parts B = Bi n + B ex . The first part of the integral Bi n is taken over the 



range \z — z\ < A, where a « A << £l being the characteristic length scale of the 
function X(z). The second part B ex is taken over the rest of the wire length A < \z — z\ < L, 
whereby the volume integral reduces to the integration over the wire length only. That is, 

//o [IxR 



Bex ~47ri R 3 d1, 

where I is the total current. The first integral can be evaluated in general with the assump- 
tion of a uniform current distribution and in the approximation a << £l, \X\/£l << 1- The 
first inequality has been implicitly introduced in fact when we set the length A, the second 
inequality just represents the fact that the wire is slightly deflected from a straight line. The 
first integration asymptotically gives the expression for the net force as follows 

{F x } in = -^0{ln(2A/a) - 3/4} + 0((X/^) 2 , A/&, (a/A) 2 ). (7) 



That is the force is proportional to the local wire curvature |IT| • The second part cannot be 
integrated in a general case and one needs to make certain assumptions about the wire form 
and thus about the function X(z) itself. In the case of an infinitely long wire with periodic 
lateral perturbations X = X cos(kz), k = 27t/£l, one gets 

{F x } ex = + 1/2 - ln(*A)} + 0(kA, (X/^) 2 , (a/A) 2 ), (8) 

where C ~ 0.577 is the Eiler's constant and we have substituted Xk 2 = - ^f . Combining 

oz 

two parts we obtain 

Fx = ~^l^ {H2/ka) ~°~ 1/4} - (9) 



Note, that the expression (gD is identical to the formula obtained in J[2| in a particular case 
of small periodic deflections \X\ « a. 

It is clear that the assumption of an infinitely long wire to calculate the force, F x , is not 
very rigorous, since in experiments the circuit is always closed. This implies that in addition 
to the magnetic force by the wire current itself, there may be a component generated by 
the currents flowing in the external circuit. This effect can be taken into consideration by 
substituting an appropriate circuit form, in other words the function X(z) in the second 
integral. 



Alternatively, this force can be obtained by means of the virtual work principle. To 
estimate how strong the effect could be, we will use a general expression for the inductance 
of a closed circuit A = ln(L s /a s ), where L s is the circuit characteristic length and a s is 
the radius of the circuit wire, see ||13|| . Then varying the circuit length L s one can calculate 
the subsequent variation of the associated magnetic energy E m = \^I 2 - Then, equalising 
the net force acting on the circuit Fn with the variation of the energy we get 

F N 5L S = ^I 2 (\n(L s /a s ) + 1)5L S . (10) 

An average force acting on a unit length is / = F^/L. Then, for the average force 
from ( |TUD we get / ~ ^ o I 2 (ln(L s /a s ) + 1)/L S . The first, logarithmic part of the net force 



represents contribution from the end effects, see (0). The second part gives the integrated 
force from the circuit as a whole. If L << L s , one can assume that the force is constant 
throughout the wire length in the first approximation. Then, this " constant" force acting on 
a wire of length L will create longitudinal stress cry = (jQ4^Q ^f^lhi see the problem 



in 



T0|| on page 93. For a copper wire with L = 30 cm, a = 0.6 mm and the circuit with 
I = 5kA, L s = 20 m and a s = 1cm the stress is amounting to cry ~ 8.4 MPa. This is only 
an estimation, and the effect needs special analysis since the geometry of external circuits is 
usually unknown. Moreover, using a symmetrical circuit with two loops from both sides of 
the wire one can completely compensate the external magnetic force. 

We leave this question for further detailed investigation and neglect the contribution from 
the external circuit in our analysis. Further we will use formula (H). 

2.2 The temperature behaviour 

The temperature behaviour with time is governed by the passing current. It can be calculated 
according to the direct Joule heating of the wire material from 

dT 

pCv-^=f/(T, (11) 

neglecting by the process of thermal conductivity. This assumption would hold for thermally 
isolated wires with uniform current distribution. In expression (|Tl~]), j is the current density, 



p and a are the density and the electrical conductivity of the metal, c v is the specific heat. 

From the above equation using the well-known inverse dependence of the conductivity on 
T 

temperature a = cr -jl one can obtain 



T = T exp( 7 T f f(t') 2 dt'), 
Jo 



where 77- = — ^ — is the characteristic temperature rise-time, jo is maximum current 
density. The function f(t) defines current rise time profile. 

3 Linear stability 

To investigate possible instabilities in the system, we first carry out a linear stability analysis 
assuming for a moment that all the parameters, such as current and wire temperature, are 
constants independent of time. Even though, in a real situation, they vary quite fast with 
time and numerical methods must be involved to solve the system, our simplified analysis 
will form a basis for qualitative interpretation and understanding of the results. 

3.1 Stability of an initially straight wire 

Consider an initially straight wire by letting X = 0. General analysis of the system (|1]) 
shows that two basic types of buckling instabilities may develop. They are due to thermal 
expansion, expressed by the force G applied at the wire ends |L(J, and the magneto-elastic 



buckling instability |14| expressed by force F x . In both cases, the instability has a threshold 



character. When either the force G or the current I exceed some critical value, new stable 
states appear. The initially stable state X = then becomes unstable and buckling occurs. 

We are looking for non-trivial stationary solutions of the linearized system ([]]), (|3|), (Q). 
They are given by the equation 

EJX IV - GX" - F x = 0, (12) 

with the boundary conditions 

X = X' = at z = and z = L. (13) 



Here, the temperature, T, and the total current, J, are assumed to be independent of time, 
and so are the functions G = -aES[T - T ) and F x = - &^X" {ln(2/ ka) - C - 1/4}. 

To obtain a solution to ( |T2| ) it is convinient to use a complete set of orthogonal functions 
{Xi} , i > 1 defined by the problem 

XT + \W = o 



Xl(0) = X;(L) = X i (0) = X l (L)=0 (14) 



with the orthogonality given by 

f 0, i^k 

(15) 

r o, z ^ *; 

This set of functions is commonly used in the theory of stability and buckling of elastic 
columns |15| . 



The eigenvalues of the boundary-value problem fll4j) satisfy the dispersion equation 

1 — cos(AjL) = sin(AjL) (16) 

So, one can see that the set of functions consists of two subsets. The first subset is defined 
by the eigenvalues given by 

cos(AiL) = 1 

sin(A;L) = 0. { ' 

That is 

Aj = 7r • (i + 1)/L, i = l,3,5.... 
In this case the associated eigenfunctions are given by 

Xi = cos(Aiz) - 1. (18) 

And the second subset is defined by 

/ -i j % 4 - (AjL) 2 
cos(AiL) = - ' -L 



sin < A < L > = ITJkl) 



4+<W3' (19) 



The associated eigenf unctions are given by 
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Xi = cos(Xiz) - 1 + 2z - — sin(Ajz). (20) 

The first eigenvalue of this subset is equal to A 2 = 8.986819/L. Further eigenvalues are 
approximately given by A,; ~ n ■ (i + 1)/L, i = 4, 6, 8 . . .. 



Now, with the help of ( |lq) one can obtain from ( |12]) a criterion when the instability first 
appears 

a(T - T )ES - ^/ 2 (ln(Aia) + 0.14) > EJ\\. (21) 

47T 



The first term on the left hand side of (|21| ) is responsible for the thermal expansion effect, 
while the second term represents magneto-elastic buckling. It should be noted that while 
we are considering both effects simultaneously, for parameters relevant to wire explosion 
experiments these terms have different orders of magnitude. 

The criterion obtained shows when the first eigenmode of ([I]) and the system as a whole 
become unstable, Ai = 2ix / L being the corresponding eigenvalue. In a similar manner, 
substituting other eigenvalues Aj for Ai, one can obtain respective criteria for higher modes 
as well. 

The set fli"4] ) being very usefull in the linear stability analysis of system ([!]) doesn't seem 
to be very relevant to study the increments of the instability since 



[ L X t X k dz ^0, i^k. 
Jo 



But, it would be interesting to obtain an estimate of the increments of the instability. For 
this purpose, we will substitute a solution to (Q) in the form X = Ai(t)Xi, i.e assuming that 
only one mode is present. Then, for the increment ji, (Ai(t) ~ exp(jit)), one gets 

7? = -4^— {S-^} (22) 
Xfdz 1 S J 

S = «(T-T )-£^{0.14 + ln(aA l )} 



or since J Xfdz ~ L 



— - \f-Yl> (23) 
p S 
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As is seen from (Pg), the increment has a maximum at Aj = X ex t defined by the equation 

AL = - £^{°- 64 + MaA^)}. (24) 

As we will see further, for the range of parameters used in the wire explosion experiments, 
the contribution from the terms due to magnetic force can be neglected during the initial 
stage and within this approximation, using explicitly J = 7ra 4 /4 and S = na 2 , 



X ex t — y — (25) 



And, since Aj ~ 7r • (i + 1)/L, 



= /MT^y l _ L 

V a it 

On the other hand, the maximal Xu m at which instability may exist is obtained from ( ^3|) 
to give 

/ a(T-Tp) 

A«jm ^ 2y 2 ( 27 ) 



with the same accuracy. From (27) , since Aj ~ it ■ (i + 1)/L, one gets for the maximal mode 
number 

ium H JaV^tLk - 1 (28) 

V a 7r 

Thus the dominant mode lies somewhere between the first unstable mode and the last one 
with the maximal increment 

7- = MaHzM. (29) 

VP a 

From the expression (pq ) one can see that higher modes are likely to become dominant 
in the instability spectrum with increase in the wire length. As a result, the wire shape can 
take a rather intricate form. It seems that the instability on higher modes has been observed 
experimentally by Graneau, [§J. At the experimental conditions he used, i.e. aluminium 
wire, L = lm, a — 0.6 mm, from (p5|) one gets X ex t — 294 or the corresponding mode 
number i ext = 93 at T = T me i t . The photographs of the wire shape after the current had 
been switched off showed clearly the appearance of at least 5 bulges on a small part of the 
wire. 
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We will refer further to the sets of parameters presented in Table I corresponding to 
different wire lengths and different currents. One should note that two sets, [B] and [E] , are 
relevant for wires and currents used in the experiments [|J and |I] respectively. 



Tabled 





material 


a 


L 


/ 


Pi 


P2 


P 3 


A 


aluminium 


0.6 mm 


0.05 m 


5 kA 


1.19 x 10 3 


6.11 


114 


B 


aluminium 


0.6 mm 


1 m 


5 kA 


1.19 x 10 3 


13.6 


0.28 


C 


aluminium 


0.6 mm 


1 m 


2 kA 


1.19 x 10 3 


2.18 


0.28 


D 


aluminium 


0.6 mm 


0.3 m 


8 kA 


1.19 x 10 3 


34.8 


0.28 


E 


copper 


0.5 mm 


1 m 


500 A 


1.84 x 10 3 


0.14 


0.24 



For all the represented sets, the criterion fl2~T|) is fulfilled well enough. Indeed, we have 
put in the Table I the corresponding values of the terms in (pl~l), designating them as Pi, 
P 2 and P3 from left to right. The estimations have been done at T = T me i t , T me i t being 
the melting temperature; T me n = 660° C for aluminium and T me i t = 1085° C for copper. 
The characteristic increment ( p9|) at the same conditions, in the case [B] for instance, is 
jmax ~ 10 5 sec -1 , which is much greater than the temperature increment 7^ — 7 x 10 2 sec -1 . 
Thus, even though the estimation has been done in the linear approximation, it is clear that 
the instability has sufficient time to develop before the wire reaches the melting point. 

In all the cases, as is seen, the major contribution is from the thermal expansion effect, 
while the influence of the magnetic-force terms can be neglected during the initial stage (this, 
of course, can be estimated directly from (p]) as well). Two terms, P\ and P2, for instance 
in case [B], become equal only at T — To = 6.6° C. The time it takes for the temperature to 
be driven through this range is At ~ 30 fisec. The characteristic time for the instability to 
develop for this temperature difference is much longer, 7 -1 x = 800 /isec. 

3.2 Stability of an initially bent wire 

Consider a wire which has an initial form Xo = J2iA®Xi, where the functions is the 
set ( |i~4"|) and A® are given weight coefficients. One needs to stress, that with the accuracy 
the system (H) was derived, \Xq\ > a must always be the case. Performing analysis similar 
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to that leading to equation (^j]) and neglecting the magnetic force for the sake of simplicity, 
for a mode k one gets 

\ 2 k (A k - A° k )EJ - A k ES(aAT) = 0. (30) 

From ( |30"D one can see that there is always a stationary state which differs from X at any 
nonzero value of the parameter aAT. Thus, those modes which have A k ^ are always 
unstable. The result is not surprising, since from the physical point of view it is obvious 
that a constantly heated wire will change its form owing to expansion. On the other hand, 
those modes which have A k = become unstable at some value of the parameter aAT, so 
that buckling instability is still possible on these modes. 

4 Qualitative stress analysis. 

Once the instability occurs, all the potential energy comprised into the compressed wire can 
be quickly released. The characteristic value of the longitudinal stress r zz which may be 
accumulated during the preconditioning is within the interval between the maximum value 
T ™ x = a(T meU - T )E ~ 10 3 MPa and the minimum value r™ n = E^\% ~ 0.3 + 30MPa 
corresponding to the onset of the instability. The estimation has been done for an aluminium 
wire of radius a = 0.6 mm and the length spanning the range lm>L> 0.05 m. 

One can see, that the maximum value itself is very high, about 10 times higher than the 
ultimate strength value. On the other hand, the minimum value can be quite small depending 
on the wire length. The average transverse stress is given by <r xx >= Q,X' - E^X'". For 
small deflections, <C 1, it is in general smaller then the longitudinal one. 

From the qualitative point of view it is evident that the maximal stress energy can be 
accumulated if the temperature rises sufficiently quickly. The temperature rise time must be 
shorter or comparable with the rise time of the instability. The most dramatic result might 
be expected if the wire is heated up to the temperature just slightly below the melting point. 
Then the current is switched off thus allowing the instability to develop without melting the 
material. On the other hand, if the current is low and, as a consequence, the temperature 
rise time is low too, then all the accumulated energy can be released at the onset of the 
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instability by the first mode which is becoming first unstable. This case corresponds to the 
lower limit of the estimated stress value r™ m . 

Thus, already now, from the linear analysis, it is obvious that many scenarios of the 
instability are possible. Dynamically, many modes can be excited simultaneously while 
the temperature T increases from To up to the melting point T me i t . Even though the first 
mode becomes first unstable it might happen that further other modes play dominant role 
creating rather complex dynamical behaviour. Any particular pattern, of course, depends on 
temporal characteristics such as the ratio between instability increment and the temperature 
rise time. 

5 Numerical results 



Now, there are several further questions. What character will the instability have on the 
nonlinear stage? How high longitudinal and transverse stresses could be obtained in this 
situation? To answer these questions we solve the nonlinear equation ([]]) numerically. 

5.1 Numerical algorithm 

The numerial method of the solution of equation ([!]) with the boundary conditions (H) is 
based on the complete set of orthogonal functions fllij) . 

By expanding any solution to ([[]) into the series of the functions {Afj}, i.e X — J2i AiXi, 
the partial differential equation turns into a system of ODEs with respect to time. They are 

E Bij^- + A~ + A~ {(3 a (l-l - a[T(t) - T ]) + ft, (0.14 + ln(A,a/L))} = (31) 



dr 2 2 



where the coefficients are f3 a = 4L 2 /a 2 , (3b = 4^gjj • ^he symmetric matrix B^ is given by 
Bij = Jq 1 XiXjdz. The length increment is I — Iq = ^ Jlj(A? — A® 2 )\ 2 . The functions X, X 
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and the variable z have been normalised on L and t has been normalised by to = J 
In the numerical simulations, the infinite series have been trancated at the maximal unstable 
mode defined by (|28D . 
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5.2 Flexural vibrations of initially straight wires 

In the case of initially straight wires the initial conditions are A(0) = 0. To excite the 
instability one needs to seed some initial noise level. To simulate noise present an additional 
fluctuation force in the equation (|1|), 5F = 6F Sf(t) Y,i ex P(~ £f/£i)Xi, has been added, with 
SF = 4 x 10 -4 newton. The function —1 < 8f(t) < 1 has been calculated by means of a 
random number generator at each step of calculations over time. This force can give rise 
to deflection of an aluminium wire with L = lm and a = 0.6 mm from a straight line by 
approximately X ~ 0.1 a. 

The current profile in the simulations has been taken in the form 
I(t) = Iq sinh(t/to)/cosh(t/to), where to is the current rise time. In all runs to = 30 //sec. 
This time is about twice as high as the skin time for aluminium wires with a = 0.6 mm, 
tskin = aV/io — 17/isec. We also kept T = 300° K and calculations stopped once the tem- 
perature had reached the melting point T = T me i t . We designated this moment by t = t e . 

We have performed simulations for different conditions: different wire lengths, different 
currents and different wire materials. All the sets used are presented in Table I. 

As the first example, let's consider a short aluminium wire carrying 5kA current, case 
[A] . For these parameters from the linear analysis one would expect a few unstable modes to 
develop, iu m = 6. The results obtained are presented in Fig.l. We have plotted in the first 
two frames the deflection A as a function of z and the spectrum < A(t)j > both taken at t — 
t e . In the last two frames, we have plotted the longitudinal t zz and transverse < t xx > stress 
components as functions of time. The spectrum has been calculated by means of averaging 
over time around t = t e , i.e. < A(t e )f >= // e e _ A A(t)fdt, where A is chosen to be greater 
than the period of nonlinear vibrations. Both the spectrum, Fig. lb and the wire shape, in 
the form being close to a simple arc, Fig. la, show that the first mode dominates during the 
nonlinear stage of the instability. Temporal behaviour of the longitudinal and transverse 
stress components demonstrates developed nonlinear flexural vibrations, Fig.lc and Fig. Id. 
The appearance of the instability is clearly seen at t ~ 900/isec in both figures. The observed 
longitudinal stress has both compressive and tensile components. The instability appeared 
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after strong compression, the maximal compressive stress being 429 MPa. The maximal 
tensile longitudinal stress is amounting to r zz ~ 244 MPa. This value is well above the 
ultimate stress value for aluminium. The observed transverse stress, as has been expected, 
is lower then the longitudinal one, < r xx >~ 55 MPa. 

Let's consider now a longer wire, case [B]. This case is relevant to the conditions of the 
Graneau's experiments ||J]. One might expect, in accordence with our qualitative analysis, 
more active modes in the instability spectrum to develop with respect to the previous case, 
him = 129. Indeed, from Fig. 2a and Fig.2b one can observe that the spectrum becomes 
very rich with the major contribution from mode n = 27 at the melting point. The tensile 
longitudinal and transverse stresses are amounting to t zz = 122 MPa and < r xx >= 92 MPa 
in this case, see Fig. 2c and Fig. 2d. Both of them are well above the ultimate strength 
of aluminium. Thus, one might expect the first wire break to appear just after buckling 
occured, at t = 1000 fisec. 

Let's investigate now the effect of variations of current. In case [C], which is a replica of 
case [B] but with lower current, two modes n = 5 and n = 7 become dominant, Fig. 3a. This 
fact indicates, if we compare two spectra in cases [B] and [C], that the main channel of the 
accumulated energy to release is going now through the modes with lower numbers. Thus, 
one would anticipate a lower accumulated stress according to our qualitative stress analysis. 
Indeed, we have observed that both tensile longitudinal and transverse stresses reduced to 
t zz = 70 MPa and < t xx >= 23 MPa, Fig.3b and Fig.3c. 

On the other hand, an increase of the current leads to higher resultant tensile stress 
value. In case [D] which is similar to set [B] but with a higher current, the longitudinal 
tensile stress is r zz ~ 233 MPa in maximum, Fig. 4. It should be noted that because of faster 
heating the preliminary compression before the instability developed is much more stronger, 
namely 851 MPa, and the instability occurred just before the melting point. It appears that 
at higher current the wire will be molten before the buckling occurs. 

Let's consider now the Nasilovski's experiment with a long copper wire, case [E]. This 
experiment has been carried out at quite a low current / = 500 A in comparison with the 
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all previously considered cases. As a result, one might expect low stress values compared to 
previous cases. Also, a few first modes must be dominant. Indeed, that is just the case, see 
Fig. 5a. The observed values of stress indeed become quite low. The longitudinal stress is 
reaching t zz ~ 27MPa in maximum, Fig.5b, while the transverse stress is just about 4MPa. 
It means that the wire can be broken only if it was heated up enough. 

5.3 Flexural vibrations of initially bent wires 

In the above, we considered somewhat idealised situation, when the undeformed wire had the 
shape of a straight rod. In reality, the wire might have had initially any form. For instance, 
in the case of horizontal positioning of the wire it might be bent by gravitational force. Once 
the initial dislocation becomes grater than the radius, the character of the instability changes 
as has been discussed. Indeed, if X ^ the initial state is unstable from the beginning. 
But, if the temperature risetime is shorter than the period of flexural vibrations then those 
modes which contribute into X = J2i A®X{ can be excited directly. Moreover, buckling on 
the other modes is still possible. To estimate the stresses developed simulations for case [B] 
have been performed but with initial conditions given by X = A\X\ with A\ = 20a. That is, 
the wire was initially shaped like the first mode of the set fllil) . The choice of the amplitude 
A\ follows from estimations for a horizontally positioned wire of this length. The wire in 
this case has a form X = — 24Ej"^ ' see ^ e problem in []TU[ page 93, q being the wire 



weight per unit length. Thus, the maximal displacement is expected to be X 0max ~ 18a. 

The results of simulations are shown in Fig.6. As before, one can see developed nonlinear 
vibrations. The regular character of the vibrations points out to the fact that only one 
single mode i = 1 is active in this case. This means that buckling on the other mode didn't 
occur. The amplitude of tensile longitudinal stress is amounting for this particular case to 
t zz ~ 50 MPa. So, a sufficiently heated wire can be broken in this well. 
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Conclusions 



When an electric current passes through a thin metal wire with clamped ends, flexural elastic 
stress waves are induced owing to the Joule heating and the electromagnetic force. The Joule 
heating leads to thermal expansion of the wire material, which is the dominant mechanism 
of the excitation of vibrations. Under realistic experimental conditions the electromagnetic 
force is of minor importance. 

Flexural vibrations in an initially straight wire may be excited as a result of the buckling 
instability. The energy accumulated in the wire during the initial stage in the form of a 
compressive stress is suddenly released. As a result, high tensile stress appears, which is 
sufficiently high to cause the fracture of the wire within 1 ms for all the cases considered. 
The number of modes of the instability that are excited depends on the current magnitude 
and wire length. Depending on these parameters it is possible to excite just a single mode, 
which has a form of an arch, or actually any number of modes. 

If a wire is slightly curved, which is a more realistic case than that of the straight wire, 
the buckling instability is still possible. The amplitude of modes, that are initially present in 
the curved wire, grows. The other modes can still be excited in a rapidly heated wire owing 
to buckling instability. The magnitude of tensile stress induced in a curved wire is clearly 
lower, but is still sufficient to induce a wire fracture on a millisecond timescale. 

Three-dimensional effects, which may be induced by the external circuit, suspensions at 
the clamped ends, imperfections of the wire cross-section or a wire material, will increase 
the magnitude of tensile stress. These effects will lead to the coupling between flexural, 
longitudinal, and torsional modes. This study is beyond the scope of the present paper. 

The experimental evidence, in general, is supportive of the mechanism described above. 
However, no direct comparison is possible, since previous experiments were of exploratory 
nature. 

The model developed favours the view that the phenomenon of the wire fragmentation 
in the solid state can be explained without resorting to controversial Ampere-Neumann 
electrodynamics. 
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Captions 

Fig. la The wire displasement X/L as a function of z/L, case [A]. 

Fig. lb The spectrum of the wire vibrations at the melting point < A} > as a function of 
mode number i, case [A]. 

Fig.lc The longitudinal stress r zz as a function of time, case [A]. 

Fig. Id The transverse stress < t xx > as a function of time at z/L = 0.9 at z/L = 0.9, 
case [A]. 

Fig.2a The wire displasement X/L as a function of z/L, case [B]. 

Fig. 2b The spectrum of the wire vibrations at the melting point < A\ > as a function of 
mode number i, case [B]. 

Fig.2c The longitudinal stress t zz as a function of time, case [B]. 

Fig.2d The transverse stress < r xx > as a function of time at z/L = 0.9, case [B]. 

Fig. 3a The spectrum of the wire vibrations at the melting point < A\ > as a function of 
mode number i, case [C]. 

Fig. 3b The longitudinal stress function of time, case [C]. 

Fig.3c The transverse stress < r xx > as a function of time at z/L — 0.9, case [C]. 

Fig.4 The longitudinal stress t zz as a function of time, case [D]. 

Fig. 5a The spectrum of the wire vibrations at the melting point < A} > as a function of 
mode number i, case [E]. 

Fig. 5b The longitudinal stress 

l~zz function of time, case [E]. 
Fig.6 Flexural vibrations of an initially bent wire, longitudinal stress as a function of 
time. Xq = A\Xi, A\ = 20 a. The wire and current parameters are relevant for case [B]. 
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